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Abstract 

We derive an expression for the free energy of the blends of block copolymers expressed 
as a functional of the density distribution of the monomer of each block. The expression 
is a generalization of the Flory-Huggins-de Gennes theory for homo polymer blends, and 
also a generalization of the Ohta-Kawasaki theory for the melts of diblock copolymers. The 
expression can be used for any blends of homopolymers and block copolymers of any topological 
structure. The expression gives a fast and stable computational method to calculate the micro 
and macro phase separation of the blends of homopolymers and block copolymers. 

1 Introduction 

Block copolymers and their blends show various interesting micro structures in equilibrium due to 
the difference in the segmental interaction between the polymer blocks pQ. The phenomena has 
been extensively studied by the self consistent field (SCF) theory |21 01 [HI IH1 IZ1 IH1 [HI [TU1 [ZH in 
which the free energy of the system is calculated by evaluating the path integral (or by solving 
the Edwards equation) for each polymers in the mean field. This method is quite useful as it can 
calculate the free energy of any blends of polymers of arbitrary topological structure. On the other 
hand, the SCF calculation is computationally demanding as it needs large memory and large CPU 
power. 

Many works have been done to reduce the CPU time and the memory needed for the SCF 
calculation. If the system has a periodic structure of known symmetry, the equilibrium structure 
(or the local equilibrium structure for the given symmetry) can be calculated rather efficiently by 
obtaining the eigen functions for the Edwards equation E] > or by using the narrow interface 
approximation and the unit cell approximation 00]. On the other hand, if the system does not 
have periodic structure (such as in the case of micellar systems) or if the symmetry of the system 
is not known, other strategies are needed to get equilibrium structures. Fraaije and his coworkers 
developed a theory for the dynamics and conducted several dynamic simulations, such as the 
dynamic behavior of the block copolymer melts under shear flow or the aqueous solution of the 
block copolymers [HlEHZj- Fredrickson et. al. proposed an efficient SCF algorithm for real space 
simulation [H], and introduced complex fields to improve the efficiency of convergence 

Though significant progress has been made in the SCF calculation, it is still very difficult to sim- 
ulate the 3 dimensional large systems of reasonable size of polymers as the required computational 
resources increase dramatically with the increase of the chain length: for 3 dimensional simulation, 
the calculation is practically limited to the system where the effective degree of polymerization is 
not so large. 
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An alternative approach is the density functional (DF) theory first introduced by Leibler |12|. 
In this approach, the free energy F of the system is expressed as a functional of the density 
distribution function of each monomer species. Clearly the computational cost of such approach is 
much less than that of the SCF calculation since the evaluation of the path integral which costs the 
major part of the cpu time and memory in the SCF calculation is not needed. For example, consider 
to simulate the block copolymer solutions: the degrees of polymerization of solvents are unity and 
ones of block copolymers are large. We cannot assume any periodicity for such a systems, thus 
we need to do the real space simulation. The computational cost of the SCF calculation increases 
dramatically as the degrees of polymerization of block copolymers increase, while the cost of the 
DF simulations is independent of the degree of polymerization. 

In the first DF theory [21 applied to AB type diblock copolymers, the density functional was 
obtained as a power series of 5(f>(r) = </>(r) — (j>, where cj)(r) is the volume fraction of A segment at 
point r, and 4> is the spatial average of 4>{r). 

0=y J dr<t>(r) (1) 

where V is the volume of the system. The free energy is expressed as ^U] 

F[<t>(r)]=F[$\ 

drdr'T(r - r')6<p(r)6<p(r') 

I / (2) 
+ — drdr'dr"T {:i) {r -r',r- r")6<p(r)6<p(r')8<t)(r") 

-I 

where F [</>] is the free energy for homogeneous state. The vertex functions T(r), and T^(r,r') 
are calculated by solving the Edwards equation or by calculating the density correlation functions. 

Leibler's method is mathematically sound, but it can be applied only for the the case of weak 
segregation, i.e., the case that the deviation from the homogeneous state is small. In order to 
apply the DF theory to the strong segregation case, Ohta and Kawasaki proposed an approximate 
expression for the free energy ^J^J. They noticed that in Leibler's theory, the term which is 
essential for giving the micro phase separation is the "long range" term. This term appears in the 
limit of small wave vector q in the Fourier transform of T(r) 

T(q) = J dre- tq r T{r) (3) 

For small wave vector q, T(q) diverges as follows: 

r(q) = 4 for « 2 ^0 (4) 
q- 

where A is a certain constant. The divergence comes from the fact that the monomer A of block 
copolymers are always connected to monomer B and therefore cannot separate from monomer B 
further than the size of the block polymer. Ohta and Kawasaki then proposed the following form 
of the free energy functional for the diblock polymer: 

F [0(r)] = F L [0(r)] + F S W)} (5) 

where Fl \4>{r)\ and F$ [<fi(r)] stands for the free energy due to the long range and short range 
interactions. The long range part is expressed as 

F L [<t>(r)} = J drdr' -Q{r - r')ty(r)<ty(r') (6) 

where Q(r) is the inverse Fourier transform of 1/q 2 . 
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where d is the dimension of the space. In the real space, Q(r — r') satisfies 



V 2 G(r-r') =5{r-r') 



(8) 



As to the short range part, Ohta and Kawasaki assumed the Cahn-Hilliard free energy: 



F s [<t>(r)] = dr -|V0(r)| 2 + /(^(r)) 



(9) 



where B is a constant and / (</>( r )) is a certain function which has double minima. 

As the Ohta-Kawasaki theory involves a few intuitive arguments, generalization of the theory 
to other systems was not straightforward. Indeed the generalization has been done only for limited 
systems such as the ABC triblock copolymers ^JE], the blends of AB diblock copolymer and C 
homopolymer |17j . and the blends of A homopolymer, B homopolymer and AB diblock copolymer 



In this paper we propose a general expression for the free energy which can be applied for any 
types of block copolymers and their blends. All parameters in the free energy are expressed by the 
microscopic parameters appearing in the SCF theory, i.e., the structure of block copolymers and 
the interaction parameters (the x parameters). This expression gives a fast computational method 
to calculate the equilibrium structure formed by the blends of block polymers. 

An approach similar to ours was already taken by Bohbot-Raviv and Wang |19|: they pro- 
posed an approximate expression for the free energy functional which can be applied for arbitrary 
structure of block polymers. We will compare our theory with their theories later. 

The present paper is constructed as follows. Firstly, we obtain the free energy functional 
for the situation that the deviation from the homogeneous state is small. This step is rather 
straightforward: the free energy is obtained by using the standard procedure of the linearized 
mean field approximation, often referred to as the random phase approximation (RPA) |201 HOj . 
Secondly we generalize this expression and seek an expression which can be applied for the situation 
that the density deviation from the homogeneous state is large. The second step is rather arbitrary, 
but we found an expression which reduces to the Flory-Huggins-de Gennes theory for homopolymer 
blends |23 an d to Ohta-Kawasaki theory for the melts of diblock copolymers. Thirdly we test this 
expression of the free energy by comparing the results of this expression with that of the SCF 
theory, and discuss the validity of the expression. 

2 Free Energy Functional for Small Density Perturbation 

2.1 The Second Order Vertex Function 

We consider a mixture of homo polymers and block copolymers which consist of several types 
of monomer units A, B, C... The chemical structure of block copolymers are characterized by 
the species of monomer units (such as A,B,C...), the number of monomers in each block, and the 
connectivity of these blocks. Polymers having the same chemical structure, i.e., polymers consisting 
of the same blocks (of same monomer units and of same length) connected in the same way are 
regarded to belong to the same polymer type. The polymer type is distinguished by the suffices 
p,q,r, . . . , and each block in the polymer is distinguished by the suffices, i,j,k,... (see Figure 0. 

Let Npi be the number of monomers belonging to the the i-th. block of polymer p. The total 
number of monomers in the polymer p is given by 



CHI 




(10) 



The block-ratio of the i-th block of polymer p is defined by 




(11) 
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We assume that all monomer units A,B,.. have the same specific volume and same bond length 
(this assumption can be removed but here it is assumed for the sake of simplicity). Let <fr P i(r) be 
the volume fraction of the monomer belonging to the block (p,i) at point r. Our objective is to 
find out the free energy expressed as a functional of {4> P i{r)} ■ 

We first consider the situation that the system is homogeneous at equilibrium, and calculate 
the change of the free energy for small density variation. In the homogeneous state, <j)pi{r) is equal 
to 

4> P i = f P i4>p- (12) 

Let 54> P i(r) = 4> P i(r) — f P i4> p be the deviation of the monomer distribution from the homogeneous 
state. The free energy functional for the system can be expressed as a quadratic form of 5tfi P i(r) 

H2E21 

F[{<t> pi (r)}} = F[{f pi t P }] 

+ \Y1 J drdr'T piM {r -r')8<P pl (r)5<p qj {r') ( 13 ) 



where F \{f P i(j) p }\ is the free energy for the homogeneous state, and T pi ^ q j(r — r 1 ) is the second 
order vertex function. In the Fourier space, eq l|13|) can be written as follows 

F[{<j> pi {q)}] = F[{f pi $ p }] 

+ 2 XI T^)d J d ^pi,gj(q)S(l)pi(-q)S<p qj (q) (14) 

P i,qj 
-I 

where S(f> p i(q) and ^ P i,qj(q) are the Fourier transform of 5<f) p i(r) and T p i^ q j(r): 

<5<M<?)= I ' dre-^S^r) (15) 



r p *,« (q) =Jdr , "'"V.,,, , (r) (16) 

The second order vertex function Y p i^ q j(q) determines the density fluctuation in the homoge- 
neous state and can be related to the density correlation function S P i^ q j(q) which is defined by 

S P i,M = (S<P pl {q)5<P qj (-q)) (17) 

where (...) stands for the ensemble average for the equilibrium state. The density correlation 
function S p i_ q j(q) is related to Y P i,qj(q) by the following equation: 

'^^^ P i,qj( ( l)Sqj : rk(q) = S pr Sik (18) 

q,j 

For ideal chains for which there is no interaction between monomer units, the density correlation 
arises from the connectivity of the polymer chain, and exists only for the monomers belonging to 
the same polymer chain. Hence S P i^ q j (q) for the ideal chain is written as: 

S^f ) (l) = ^ Pq h P ,i j (Q) (19) 

where h p ^j(q) is the correlation function for a single chain. Let r p f be the position of the s-th 
monomer in the i-th block of the polymer p, then h p ^j (q) is given by 

1 f N pi f N p] 

hp,a (g) = jfJ Q ds J o ds ' ( ex P h ' (r S " r S )J ) ( 2 °) 
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For Gaussian chain, the average can be calculated easily to give 

+ U (i = j) 

W«) = < Jf.f . „ „ ( 21 ) 



where 



6 is the effective bond length and 



Zpi = \N p f pl b 2 q 2 (22) 

l 2 Ptij = M p , l3 b 2 (23) 

Here M p .y is the number of monomers included in the shortest path connecting the i-th block and 
the j-th block of the p-th polymer, i.e., the chemical distance between the i-j blocks. 

Thus for ideal chains, the correlation function can be calculated by eqs 1)19(1 and ((21(1 . and the 
vertex function can be calculated by eq i(18|) . This gives the following expression for the vertex 
function of the ideal chain T^f^ (q): 

r £f(<?) = ¥hp,M ( 24 ) 
<P P 

where g P ,ij(q) is the inverse of the matrix h p ^j(q): 

^2 h p ,ij {q)9p,jk (<?) = $ik (25) 

3 

The effect of the interaction among the monomer units can be taken into account by adding 
the following interaction energy term 

Fi n t [{4>pi(r)}} = 2 X] / dr X P Lqj^pi{r)<t) q j{r) (26) 

Pi A3 

where Xpi,qj ^ s the X parameter for the interaction between the monomers in the block (p, i) and 
those in the block (q, j). Eqs 124(1 and 1(26(1 give the following vertex function: 

Tpi A i (q) = -^-g Pyi j (q) + X P i,q3 ( 27 ) 
<Pp 

2.2 Asymptotic Behavior of the Vertex Function for Large and Small 
Wave Vector 

Though T p i^ q j(q) can be obtained by the formula given above, it is desirable to have an analytical 
expression for T pitq j{q). As it was first noticed by Leibler J2|> the characteristics of ^ P i, q j{q) is that 
it diverges as l/q 2 for small q. Ohta and Kawasaki utilized this fact in deriving their approximate 
density functional for the free energy of diblock copolymers. The divergence of T p i^ q j(q) for small q 
reflects the topological structure of the block copolymer. The vertex function r pii9J (q) also diverges 
for large q region in such a way as T p i tq j(q) oc q 2 . Our strategy of constructing the approximate 
free energy functional is to get an expression which correctly describe the asymptotic behavior of 
T p i tq j(q) for large q and small q region. Therefore we first study the behavior of g P ,ij(q) for large 
q and small q region. 

For large q, the behavior of g P ,ij(q) is easily seen. In the limit of q 2 — > oo, the correlation 
function h Pt ij{q) is totally determined by the local structure of the polymer. From eq ((21(1 . it 
follows 

W<7) = ^|^ + ---, <? 2 -°° (28) 



5 



Thus 



is given by 



Sjjb 2 
12/, 



p< 



The behavior of g Pt ij(q) for small q region is less obvious. In the limit of q 2 
be expanded into the power series of q 2 . 



hp,ij(,q) — Npfpifp 



q 2 ^0 



where the expansion coefficient H Pt ij is given by 

1 



H„ 



—N 2 f 3 b 2 

Npjpifpj 



^Np{f pi + f P] )b 2 + \ll l3 



(29) 

0, eq (|2U can 

(30) 

(31) 



Eq i|5U|) indicates that det(h Pii j (q)) becomes zero at q 2 = (since det(N p f pi fpj) = 0). Therefore 
the matrix equation H25|) becomes singular for q 2 — ► 0. It is shown in the Appendix 1X1 that the 
solution of eq (|25|) for small q can be written as 



9p,ij(q) 



A 



:>■<:] 



(32) 



where 



( H P 1 ) lk Ukfpl (Hp *) 



.4 



PM 



(Vp- 1 ), 



k! 



Yfpk{.H p 1 ) kl fpi 
kl 

where (H p ~~ 1 ) i . is the ij component of the inverse matrix of (Hp)-: 

Z2 H P,ij (Hp 1 )^ = S ik 



(33) 



(34) 



2.3 Approximate Expression for the Vertex Function 

Having seen the asymptotic behavior of g P ij (q) , we now seek an approximate expression for g p y (q) 
which can be used in the entire region of q. Considering the asymptotic behavior in the two limits 
of q 2 — > and q 2 — > oo, we use the following expression 



g P ,ij(q) = 



A 



p-i-j 



Sab 2 , 
_y — q 2 



12/ 



(35) 



pi 



where the constant C Pt ij is chosen so that eq (|35|) gives a good approximation in the intermediate 
region. We use the following procedure to determine C p ^j. 



For i = j, g P ,ij(q) has a minimum at q* pi = \J 12 fpiA Pt n/b 2 . We thus determined C Pi u so that 
the minimum value agrees with the exact value. This gives 



Cpsi — (h p 1 (q*pi)) ii 



A 



Pi 11 



*2 

q pi 



b 2 ,2 

— q P i 



12/ 



(36) 



where (h p 1 (q* P i))- stands for the ij component of the inverse matrix h p .ij{q) at q = q* pi . 

For i ^ j, g P ,ij(q) is a monotonically decreasing function of q 2 and approaches to a constant 
value Cp^j for q 2 — > oo. We thus determined C Pj ij by (/i~ 1 (co))ij. This value is given by (see 
Appendix EJ 

^Npfpifpj ^ ^ (37) 
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Cp,ij — 



G 



From eqs l|27[) and (|35[1 . the second order vertex function is obtained as 



Sijb 2 2 
-q 



Xpi,i 



This gives the following free energy in the real space representation: 
F[{^(r)}] = F[{/ pl p }] 

+E/^%^^w^w 

+ E Jdr*f>5cj> pi (r)64, qj (r) 



(38) 



(39) 



pi, 13 



This free energy functional is valid for |<S0p»(r)| <C 1. 

2.4 Test of the Approximate Free Energy for Small Perturbation 

We now test the accuracy of eq (j3"§)) by calculating the density correlation function of block copoly- 
mers. The density correlation function S P i tq j(q) can be calculated from Y p i^ q j(q) by solving eq (JTSJ. 
For polymer melts, an additional constraint called the incompressible condition, is usually imposed. 
This condition is written as 



E = or ' E Ml) = for 9^0 



p,2 p,l 

This condition is equivalent to assume the following interaction energy 



X 



E <M<2)<M-<?) 



(40) 



(41) 



and taking the limit of \ ~ * °°- This gives the following density correlation function: 



S phqj (q) = (r- 1 ) piM (q)- 



E( r_1 W*te) 



r.k 



Er r 



qj,rk 



(<?) 



r, k 



E ( r ~V,s*(g) 



r.s.k.l 



(42) 



The second term comes from the incompressible condition. 

Figure [3 shows the scattering functions for AB diblock polymer melts for xab — 0. Here S{q) 
is defined by 

S(q) =S AA (q) + S BB (q)-2S AB (q) (43) 

The dashed line denote the exact scattering function calculated by eqs l|18|) . (|27fl and (|43() . The 
solid line denote the result of the approximate scattering function calculated by eqs l|18(l , (|38J) and 
Ij43(l . It is seen that the agreement is quite good. 

Figure shows the scattering functions for symmetric ABA triblock polymer melts. The agree- 
ment is again quite good. 
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3 General Free Energy Functional 



3.1 Free Energy Variation for Phase Separated State 

Eq (|39|) can be used only for the homogeneous state where |<50 P i(r)| <C 1. We now seek a general 
expression for the free energy which can be used for the phase separated state. The expression 
has to reduce to eq (|39(l for the case of small deviation from the homogeneous state. Clearly 
such generalized expression is not unique. The correct way to handle such systems is to evaluate 
the correlations under inhomogeneous density profile, but it is quite difficult to evaluate them 
analytically (in most cases, it must be done numerically). In the following, we propose a plausible 
form of the generalized expression based on physical argument instead, and test its validity for 
some typical cases. 

First we note that, for phase separated state, eq l|39|l is not valid even if the density deviation 
54> p i(r) from the equilibrium state is small: for macroscopically phase separated system, <f) p in 
eq (I5§t should not be the average of <p p (r) for the entire volume; it should be the average of 
4> P {v) in the local region of the phase separated state where the point r and r' are located. This 
consideration suggests the following replacement 

4> p -> lMr)Mr')} 1/2 (44) 

However such replacement does not give an analytically tractable form for the free energy. As an 
alternative, we used the following replacement: 



~ <M* , )'MO j 1/2 

fpifpj 

This gives the following free energy functional for small variation of <fi p i(r). 

S^F[{d> pi (r)}} = 

Y I drdr' -^f-f-A p tf g( r - r>) 

v f P \VS^(r)\ 2 
+ J2jdr*f>6<t> pi (r)64> qj (r) 



(45) 



(46) 



3.2 Free Energy for the General Case 

The right hand side of ea. l|46ll represents the second order functional variation of the functional 
F [{4>pi( r )}) with respect to <p p i(r). The equation is then regarded as a second order functional 
differential equation. From the equation, the functional F [{4> p i (r)}] is uniquely determined. The 
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result is 



F[{<t>pi(r)} 



^2 / drdr' 2^J fpif P jA Ptij g(r - r')yj <j> 

pi { r )4>pj 



pi 



+ Y1 j dr -\ fiJ'mC,,^ sj <j> pl {r)<j) P j{r) 



(47) 




It is straightforward to check that the second order functional variation of eq (|46|l with respect 
to 6<fcpi(r) leads eq (|39H . This free energy functional does not depend on (j> p and can be used for 
strong segregation. 

Of course there ane many other possible free energy functional which reduces to eq (|39[1 , but 
eq (|4*7)l has several advantages. 

1. For homopolymer blends, eq l|47|) gives the Flory-Huggins-de Gennes type |22II21| (or Lifshitz 
type |23j) free energy. This is shown in section |4~T1 

2. For the ordered phase of diblock copolymers, eq (|47|l gives the results equivalent to that of 
the Ohta-Kawasaki theory |13M14| . This will be demonstrated in section l4~21 

3. The minimization of eq l|47|) with respect to 4> P i{r) usually requires numerical calculation 
for which eq (|47|l has an advantage. We found that by using the variable ?p P i{r) = J 4> P i(r) 
rather than 4> P i(r) itself, we can improve the numerical stability significantly. The detail of 
the numerical procedure is described in Appendix iBl 

4 Test of the Free Energy 

In this section we apply the free energy expression l|7t)|l to special cases, and discuss its accuracy. 
4.1 Homopolymer Blends 

For homopolymer, A p becomes identically equal to zero since eq (I33|) gives 




= 



(48) 



and there is no long range interaction. 
C v can be calculated as follows. 



c p = (MO)) 



i 



i 



(49) 



p 
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From eq (I47|) . the free energy functional for the system can be described as 

F[{Mr)}] = 

|V0 p (r)| 2 



+ E/ dr^Mr)Mr) 



(50) 



Thus our free energy functional reduces to the Flory-Huggins-de Gennes type free energy functional 
|2~T] (strictly speaking, the free energy given in ref [21] has the factor 1/36 in the second term of 
eq (|50ll . This is because in ref the factor is determined from the behavior of h(q) for small q 
region, while we determined it using the behavior for large q region. The factor 1/24 agree with 
that given by Lifshitz et. al. |23p. 

4.2 Diblock Copolymer Melts 

For diblock copolymer melts, Aij and CVj are given by 



A 



9 



" N 2 b 2 f 2 (l-f) 2 
1 



(I-/) 2 
-/(I-/) 

*(/) -1/4 ' 
-1/4 s(l-f) 



-/(I-/) 
/ 2 



(51) 
(52) 



10 Nf(l - f) 

where 6^ = bs = b, Ja = 1 — fs = f and s(/) is the function which is determined from eqs (|36|l . 

nm 

From eq (|39(l the free energy functional for the weak segregation limit will be 



F [#r)] 



drdr' 



2N 2 b 2 f 2 (l - f) 2 
dr x5(j) 2 {r) 
b 2 



Q{r - r')5(j)(r)5<f>{r') 



(53) 



dr ■ 



24/(1 - /) 
W[5cf>(r)} 



\V6<j)(r) 



where we set 0(r) = M*) = 1—4>b(t), F [4>] = and* = X AB-[s(f)+s(l-f) + l/2]/2Nf(l-f). 
W [5<fr(r)] is the contribution of the higher order terms for S(f>(r). In this case the free energy 
functional reduces to the Ohta-Kawasaki type free energy functional |13| . 

For strong segregation, we get the following free energy functional from eq (|47|l . 



F [#r)] = 

drdr' 

dr ■ 



18 



N 2 b 2 f(l - f) 



1 



Q(r — r')T(r)T(r') 



N 



S ^ cj){r) In <f>(r) + ^ ^ (1 - </>(r)) In (1 - 0(r)) 



dr 



1-/ 
1 



ivv7(i - 7) 

6 2 



24^(r) (1 - 0(r)) 
drxH r ) (1 - </>(»")) 



/ 

^(r)(l-0(r)) 



|V0(r)| 



(54) 
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where r(r) = — f)</>(r) — \f f (1 — <Pi r )) and X = Xab- The interfacial energy in eq (|51Jl (the 
fourth term) does not depend on /. This is consistent with the improved form of the Ohta-Kawasaki 
theory |14j . 

Figure 0] shows the phase diagram for AB diblock copolymer melts. The solid line is the results 
of the DF theory and the dashed line is the result of Matsen Bates [TTj which is based on the SCF 
theory. The critical point predicted by our DF is \N — 10.553778, while that predicted by the SCF 
theory is x-A = 10.495. The results of the DF and the SCF do not agree quantitatively but they 
agree well qualitatively. Especially notice that the DF gives a stable double gyroid phase between 
the lamellar and the hexagonal cylinder phase for small j^JV. It is difficult to get such a result by 
the previous density functional theories (the Leibler theory or the Ohta-Kawasaki theory). 

Figures [3 shows the equilibrium periods for AB diblock copolymer melts plotted against \N ■ 
Our DF theory does not agree well with the SCF theory, but it reproduces the scaling feature 
correctly: D oc N 1 / 2 for weak segregation and D oc A 3 / 2 for strong segregation. The numer- 
ical disagreement is the nature of the Ohta-Kawasaki type approximation and arises from the 
approximation for the vertex functions (eq J3HJ)- 

Figures El shows the equilibrium density profile of AB diblock copolymer melts. The solid line 
is the results of the DF theory and the dashed line is the result of the SCF theory jTTj. The result 
of the DF theory deviates significantly from that of SCF theory for weak segregation, while they 
become close to each other for strong segregation. 

4.3 A,B Homopolymer / AB Diblock Copolymer Blends 

We applied eq (|76|) to A,B homopolymer / AB diblock copolymer blends. The free energy model 
for A,B homopolymer / AB diblock copolymer melts was proposed and studied by Kawakatsu |18|. 
Our free energy functional has similar form to the Kawakatsu's theory. 

Figurc0is equilibrium structures for an A,B homopolymer / AB diblock copolymer blend. We 
set ipA — 4>b — 0.475, 4>ab = 0.05, N A = N B = N AB = 40Jab,a = Jab.b = 0.5, XAB = 1. The 
result of the SCF theory was obtained by using the SUSHI engine in OCTA system - Again the 
agreement between the DF theory and the SCF theory is not perfect, but reasonable considering 
that no adjustable parameter is used in these comparison. 

4.4 AB diblock coplymer / C homopolymer blends 

As the last example, we show the results of AB diblock copolymer / C homopolymer blends 
[T7] . These blends cause the macrophase separation as well as the microphase separation and show 
various interesting structures. Figures and |5] are equilibrium structures for AB diblock copolymer 
/ C homopolymer blends (the parameters are Nab = 10, Nc = 20, Jab, A = !ab,b = 0.5, 4>ab = 
0.2,0c = 0.8, xab = 1-2, xbc = hXcA = 0.5 and N AB = 10, N c = 20, f AB ,A = 0.35 Jab,b = 
0.65,(^45 = 0.3,^(7 = 0.7, xab = 1-75, Xbc — 1,Xca = 0.5, respectively). The simulation has 
been done for the three dimensional periodic system of 64 x 64 x 64 lattice points (the system size 
is 406 x 406 x 406). The structure of the AB diblock copolymer in Figure [S] is considered to be the 
"onion structure" observed by the experiments It is also noted that each simulation starts 

from the homogeneous state, and gives the results shown there after about 4 hours on the 3.0GHz 
Xeon workstation. It is quite difficult to get these results with such short computational time by 
the SCF simulation. 

4.5 Comparison with the other density functional theory 

In this section, we compare our theory with the density functional theory proposed by Bohbot- 
Raviv and Wang for block copolymer melts (T^|. They used the Flory-Huggins free energy for the 
disconnected blocks as the reference and took into account of the effect of chain correlation by the 
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second order vertex function with the RPA: their free energy functional is expressed as 



F [{Mr)}} =F Ief [{Mr)}} + AF [{A&(g)}] 
F rct [{Mr)}}=J2 I dr^-Mr)^Mr) 



(55) 



(56) 




(57) 



-IEj^J^^am-v) 



where 



A&(r) = Mr) ~ fi 



(58) 



Eq (|55() is exactly same as eq (j5J) except for higher order terms for A<^(r) (it can be easily shown 
by setting cf> — 1 and expanding eq H55fl by series of A^r)). 

Though their approach is simple and straightforward their theory cannot overcome one well- 
known difficulty. Ohta and Kawasaki noticed that the block ratio dependence of the interfacial 
tension calculated from the RPA is inconsistent with other theories [T2|. They then proposed to 
replace /, the block ratio, by the local density </>(r). Without this correction, the free energy 
eq (15511 cannot give correct block ratio dependence. The replacement eq (|45|l is essentially same 
as the replacement by Ohta and Kawasaki and therefore our free energy ()76|) does not have such 
difficulty. 

5 Conclusion 

We have given an expression of the free energy of blends of block copolymers expressed as a 
functional of the density distribution of monomers in each block. All parameters in the expression 
are determined by the polymer structure parameter (such as the degree of polymerization, the 
branching structure, and block ratio etc.) and the x parameter. For homopolymer blends, the 
expression is shown to reduce to the Flory-Huggins-de Gennes theory. For diblock copolymer 
melts, the expression gives results equivalent to the Ohta-Kawasaki theory. Thus the expression 
is expected to work in more general case. We have applied the expression for several cases, and 
found reasonable agreement with the SCF calculation. Of course, these examples do not guarantee 
the validity of our approximation in the general case, and more work is clearly needed to test the 
validity of our expression. 

The present theory is most useful by combining it to the SCF calculation. It uses exactly the 
same parameter set as the SCF theory (this is in contrast to the other mesoscopic approaches such 
as the dissipative particle dynamics (DPD) |27j). Therefore the theory can be used to get the 
initial structure for the SCF calculation, or it can be used to get an overview of the phase diagram 
of block copolymers. These applications will be pursued in future. 
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Appendix 



A Calculation of A p ^ and C v ^j 
A.l Calculation of A Pji j 

In this section we derive eq <|33ll . At the limit of q 2 — * 0, h Pt ij(q) can be expanded as follows. 

h p ,ij(q) = N p f pi f pj - Hp^q 2 + ■■■ (59) 

where Hp.ij is given by eq (|3ip. As we are considering one ideal chain, we omit the subscript "p" 
and describe eq (|59|) in a matrix form. 

h(q) = Nff Hq 2 + -- (60) 

where /' is the transposed vector of /. 

We obtain g(q), the inverse matrix of h(q), in the following form: 

g(g)=ii+B + - (61) 



From eq (|2*5|l . 



r/U/) • = [A±+B + --.). {Nff - Hq 2 + ■ ■ •) = E (l,2i 



and we get the following set of equations. 

A ■ f = (63) 
A H + NB ■ ff = E (64) 

where E is unit matrix. Notice that the matrix h(q) is symmetric and therefore the matrices A 
and B are also symmetric. Multiplying H 1 ■ ff to eq Q64JI from the left side and using the 
relation /* ■ A = (A ■ /)* = which follows from eq l|63|l . we get 

NB ■ ff ■ H 1 ■ ff = H 1 • //* (65) 

Since /* • H 1 ■ f is a scalar, eq <|55|) gives 

J J 1 £ £t 

NB ■ ff = £t u— i r (66) 



From eqs ijB^j) and l|6"6"jl we have 



^ = -H- 1 + ^.H-^f 1 (67) 



This gives eq l|33|) . 



A. 2 Calculation of C v ^ 

In this section we calculate the asymptotic value of g P ,ij(q) at q 2 — > oo for i ^ j. From eq (|21|l . 
^p.ij (?) can be described as follows in matrix representation. 

= G ^2 + K ^4 + ■ ■ ■ (6 8 ) 
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where 



G = 



12 /A. 
b 2 



72 

K -) 36 



M>4 rJ, y 

li -l A 12 



(* = J) 



g P ,ij(q) can be expressed as power series of q 2 . 

g (q) = J B g 2 + C + 

The equations to determine coefficients are 

BG = E 
BK +CG=0 



Thus we get 



C = BKG 1 = -G^KG 1 



(i = 3) 



2f N 



1 



WlN CM A 4-0) 



This gives eq (J37J. 

B Numerical Scheme 

The minimization of eq H47J1 with respect to 4> P i{r) usually requires the numerical calculation, 
wc briefly describe the numerical scheme we employed. 

For numerical calculation, it is convenient to introduce a new order parameter rp p i(r). 



i/j pi (r) = \J(t> P i(r) 
For this order parameter, the free energy functional is written as 
F[{i> pi (r)}} = 

^2 J drdr' 2^/J~J~A p ^g(r - r')il)pi(r)if) pj (r') 

+ dr2f pl C Piil ipl i {r)\nip pi {r) 

pi 

+ j ' dr2 \/ fpifpj C P,i3i>Pi( r ) X l ; P3( r ) 

+ E/ drj\V^(r)\ 2 



^<(r)^,(r) 



Pi, 13 
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This expression is advantageous for numerical calculation as it does not have the numerical insta- 
bility associated with the logarithmic term of eq <|47|l . The fourth term in eq (|76|l is in agreement 
with that proposed by Lifshitz |241 123j . 

The equilibrium structure is obtained by minimizing the free energy l|76|l under the following 
constraints 

J2M r ) = 1 ( ?? ) 

pi 

dr4> 2 pi (r)=Vf pi ^ p (78) 

The first equation represents the incompressible condition, and the second equation represents the 
mass conservation for each block. To take into account of these constraints we use the Lagrangian 
multiplier method, and added the following terms to eq 



constraint 



[{Mr)}] = E/ dr \ ^ + K ( r )l i^U r ) ~ fpih) (79) 



where X p i and n(r) are the Lagrangian multipliers. 

We used simple relaxation method to get equilibrium structures. 

^ +1) (r) = ^ ) (r)-^ ) (r) (80) 

where the superscript (n) means the number of the relaxation step, a; is a sufficiently small constant 
and fJi P i(r) is the chemical potential defined as follows. 

^ ^ _SF[{lp pi (r)}} t ^constraint [{lppi{r)}\ 



E< 



4 V fpif P ] A p,ij dr'g(r -r')ip p j(r') 



-^ 2 M r ) 



+ 2fpiCp t u [2ippi(r) lnip pi (r) + ip pi (r)} 

+ E A ^f fp*fp3 C P;^pj( r ) (81) 
3 

+ ^2 2 Xpi,gjM r ) 1 Pd( r ) 

13 

+ {\ P i + K(r)]ip pl (r) 

The terms including Q(i — r') or V 2 are calculated by using the fast Fourier transform (the FFTW 
library is used), and other terms are calculated in the real space. It is noted that eq l|81|) 
has no singularity at ip P i(r) = while 5F [{<fi P i(r)}} /8<f)pi(r), the chemical potential for <j) p i(r), is 
singular at (j> P i(r) = 0. Thus using the order parameter ip P i{r) rather that (f> P i(r) improves the 
numerical stability significantly. 

The numerical method we employed in this work is not so efficient We employed this method 
for simplicity, but the algorithm need to be improved (for example, the ICCG method will improve 
the convergence). 
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Figure n Chemical distance M Pt ij of block copolymers 

Figure Scattering functions of AB diblock copolymer melts. The \ parameter Xab is set to 
zero, and j a stands for the block ratio. 

Figure Scattering functions of symmetric ABA triblock polymer melts. The x parameter xab 
is set to zero, and /a stands for the block ratio defined in the figure. 

Figure 0] Phase diagram for AB diblock copolymer melts. The solid lines are the results of the 
DF simulation and the dashed lines are the results of the SCF simulation ^J. (B: BCC sphere, 
H: Hexagonal cylinder, G: Double gyroid, L: Lamellar, D: Disordered) 

Figure Equilibrium periods for AB diblock copolymer melts. Here D is the equilibrium periods 
for lamellar structure. The solid lines are the results of the DF simulation and the dashed lines 
are the results of the SCF simulation The doted line is calculated by the DF for the weak 
segregation limit. 

Figure IHI Equilibrium structures for AB diblock copolymer melts. The solid line is the result of 
the DF simulation, the doted line is the analytic solution for weak segregation limit by the DF and 
the dashed line is the results of the SCF simulation JT] . 

Figure An equilibrium structures for an A,B homopolymer / AB diblock copolymer blend 

Figure|SJ An equilibrium structures for an AB diblock coplymer / C homoploymer blend. The black 
and gray surfaces are the isodensity surface (4> P i(r) = 0.5) for the A and B segment, respectively. 

Figure An equilibrium structures for an AB diblock copolymer / C homopolymer blend. The 
black and gray surfaces are the isodensity surface (<f> P i(r) = 0.5) for the A and B segment, respec- 
tively. 
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